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Abstract 


We study the decays of J//1» and 1/(3686) to the final states (1385)99:(1385)? and E? E? based on a single baryon 
tag method using data samples of (1310.6 + 7.0) x 10° J/w and (447.9 + 2.9) x 109 (3686) events collected with 
the BESIII detector at the BEPCII collider. The decays to X(1385)99:(1385)? are observed for the first time. The 
measured branching fractions of J/4» and 1/(3686) to =°=° аге in good agreement with, and much more precise than, 
the previously published results. The angular parameters for these decays are also measured for the first time. The 
measured angular decay parameter for 7/1 — X:(1385)93:(1385)9, а = —0.64 + 0.03 + 0.10, is found to be negative, 
different to the other decay processes in this measurement. In addition, the “12% rule" and isospin symmetry in the 
decays of J/4 and (3686) to == and X(1385)3:(1385) are tested. 


Keywords: charmonium, branching fraction, angular distribution 
PACS: 12.38.Qk, 13.25.Gv 


1. Introduction 


The decays of the charmonium resonances 7/4 and 
(3686) [in the following, y denotes both charmonium 
states J/w and (3686)] into baryon anti-baryon pairs 
(BB)in e*e- annihilation have been extensively stud- 
ied as a favorable test of perturbative quantum chromo- 
dynamics (QCD) [1]. These decays are assumed to pro- 
ceed via the annihilation of the constituent cc pair into 
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three gluons or a virtual photon. 


It is interesting that the 1/(3686) decay to a specific 
final state is strongly suppressed relative to the same 
final state іп 7/0 decay according to the annihilation 
decay of heavy quarkonium. The ratio of branching 
fractions for 7 decaying into the same final states is 
predicted from factorization [2] to be Se Rd 


1296, where X denotes any exclusive hadronic decay 
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mode or ће (^/-(£ = е, p) final state. This expec- 
tation is usually called the “12% rule". This rule was 
first observed to be violated in the decay of ~ into the 
final state рт. A broad variety of reviews of the rele- 
vant theoretical and experimental results [3] conclude 
that the current theoretical explanations are unsatisfac- 
tory. Although the branching fractions for ~ decays 
into baryon pairs have been measured extensively [4], 
uncertainties are still large for many decays; e.g. the 
world average values of the branching fractions for J/w 
and (3686) — =°=° are (1.20 + 0.24) x 107? and 
(2.07 + 0.23) x 1074 [4], respectively. In particular, 
w — X(1385)93:(1385)? has not yet been observed. 

By hadron helicity conservation, the angular distri- 
bution of the process ete~ — wv — BB is expressed 
as 


dN 

dcos 0 
where 0 is the angle between the baryon and the beam 
directions in the етет center-of-mass (CM) system and 
о is a constant, which has widely been investigated 
in theory and experiment [5]. Theoretically, the value 
of а is discussed in the framework of many models, 
such as quark mass effects [6], or electromagnetic ef- 
fects [7], which generally predict 0 < o < 1. BES 
measured the angular distribution of J/y — X9»? 
and obtained a negative a with poor precision [8]. 
Н. Chen et al. [9] explained that the angular distribu- 
tion for у — BB could be negative when rescatter- 
ing effects of baryon and anti-baryon in heavy quarko- 
nium decays are taken into consideration. Thus, ex- 
perimental measurements of o are helpful to test the 
helicity conservation rule and the validity of the var- 
ious theoretical approaches. In previous experiments, 
the angular distributions for charmonium decays to 
baryon pairs, such as Y% — pp, AA, X030, £7 €*, and 
3:(1385)*3(1385)7 [10, 11, 12, 13], were measured. 
However, angular distributions for the decays v. — 
3:(1385)03:(1385)9 and E? E? have not yet been mea- 
sured. 

In this Letter, we report the most precise measure- 
ments of the branching fractions and angular distribu- 
tions for Y — X:(1385)93:(1385)? and E? E? based on 
the data samples of (1310.6 + 7.0) x 109 J/w [14] and 
(447.9 + 2.9) x 106 w(3686) [15, 16] events collected 
with the BESIII detector at BEPCII. 


x 1+acos? 0, (1) 


2. BESIII detector and Monte Carlo simulation 


BEPCII is a double-ring ete~ collider that has 
reached a peak luminosity of 10?? cm~?s~! at a CM 
energy of 3.773 GeV. The cylindrical core of the BESIII 


detector consists of a helium-based main drift chamber 
(MDC), a plastic scintillator time-of-flight (TOF) sys- 
tem, and a CsI(TI) electromagnetic calorimeter (EMC), 
which are all enclosed in a superconducting solenoidal 
magnet with a field strength of 1.0 T for the (3686) 
data and J/~ data taken in 2009, and 0.9 T for the J/¢ 
data taken in 2012. The solenoid is supported by an 
octagonal flux-return yoke with resistive plate counter 
modules interleaved with steel as muon identifier. The 
acceptance for charged particles and photons is 93% of 
the 47 stereo angle, and the charged-particle momen- 
tum resolution at 1 GeV/c is 0.5%. The photon energy 
resolution is 2.5% (5%) at 1.0 GeV in the barrel region 
(end caps regions). More details about the experimental 
apparatus can be found in Ref. [17]. 

The response of the BESIII detector is modeled 
with Monte Carlo (MC) simulations using a framework 
based on GEANT4 [18]. The production of Y resonances 
is simulated with the KKMC generator [19], the subse- 
quent decays are processed via EVTGEN [20] accord- 
ing to the measured branching fractions provided by 
the Particle Data Group (PDG) [4], and the remaining 
unmeasured decay modes are generated with LUND- 
CHARM [21]. To determine the detection efficiencies 
for) — X(1385)93:(1385)? and =°=°, one million 
MC events are generated for each mode taking into ac- 
count for the angular distribution with o; value measured 
in this analysis. The decays of the baryons X(1385)9, 
=0 and A in the signal channels are simulated exclu- 
sively, taking into account the angular distributions via 
EVTGEN [20], while the anti-baryons are set to decay 
inclusively. 


3. Event selection 


The selection of y  —  X(1385)93(1385)? 
and =°=° events via a full reconstruction of both 
3:(1385)9/Z? and 3(1385)9/5 suffers from low re- 
construction efficiency and large systematic uncertainty. 

To achieve higher efficiency and reduce the system- 
atic uncertainty, a single baryon ©(1385)°/=° tag tech- 
nique is employed, without including the anti-baryon 
mode tag due to the imperfection of the simulation re- 
lated to the effect of annihilation for anti-proton. The 
3:(1385)9 /Z? is reconstructed in its decay to m? A with 
the subsequent decays A — рт and п? — yy. The 
charged tracks are required to be reconstructed in the 
MDC with good helix fits and within the angular cover- 
age of the MDC (| соѕ0| < 0.93, where 6 is the polar 
angle with respect to the е+ beam direction). Informa- 
tion from the specific energy loss measured in the MDC 


ppl 
бе fan 


am 
e 


A 
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(dE/dx) and from the TOF are combined to form par- 
ticle identification (PID) confidence levels for the hy- 
potheses of a pion, kaon, and proton. Each track is as- 
signed to the particle type with the highest confidence 
level. At least one negatively charged pion and one pro- 
ton are required. Photons are reconstructed from iso- 
lated showers in the EMC. The energy deposited in the 
nearby TOF counter is included to improve the recon- 
struction efficiency and energy resolution. Photon ener- 
gies are required to be greater than 25 MeV in the EMC 
barrel region (| cos0| < 0.8) or greater than 50 MeV in 
the EMC end cap (0.86 < |cos6| < 0.92). The show- 
ers in the angular range between the barrel and the end 
cap are poorly reconstructed and are excluded from the 
analysis. Furthermore, the EMC timing of the photon 
candidate must be in coincidence with collision events, 
0 < t € 700 ns, to suppress electronic noise and energy 
deposits unrelated to the collision events. At least two 
good photon candidates are required. 

In order to reconstruct the 7° candidates, a one- 
constraint (1C) kinematic fit is employed for all yy 
combinations, constraining the invariant mass of two 
photons to the п? nominal mass, combined with the re- 
quirement of |AE|/ Ро < 0.95, where AE is the en- 
ergy difference between the two photons and Po is the 
1? momentum, and the x75 < 20 to suppress non-7? 
backgrounds. 

To reconstruct the А candidates, a vertex fit is ap- 
plied to all рт combinations; the ones characterized 
by x? < 500 are kept for further analysis. The рт in- 
variant mass is required to be within 5 MeV/c? of the 
nominal A mass, determined by optimizing the figure 
of merit FOM = JE where S is the number of sig- 
nal events and В is the number of background events 
based on the MC simulation. To further suppress the 
background, the decay length of A is required to be 
larger than zero. The У(1385)0 /=0 candidates are re- 
constructed with A and л? candidates by minimizing 
the variable | Mz;o, — My,1385)0/z0|, where Mron is 
the invariant mass of the т°А pair, and Ms(1385)o /=0 is 
the nominal mass of X(1385)0/Z. 

The anti-baryon candidate 4(1385)°/=° is inferred 
by the mass recoiling against the selected 7° A system, 


Mis! = \/ (Ecm — Exon)? — Poa, (2) 


where Eo, and p;o, аге the energy and momentum of 
the selected 4? A system, and Ecm is CM energy. Fig- 
ure | shows the scatter plot of Мол versus Mises 
Clear accumulations of events corresponding to the sig- 
nals of > 3:(1385)93:(1385)? and =°=° decays are 
observed. The distributions of Мо д with the additional 


requirement of the M¢¢°!! within +80 MeV/c? around 
Mys(1385)0 or £50 MeV/c? around Mzo are shown in 
Fig. 2. Clear ©(1385)°/=° signals are observed. 


(GeV/c’) 


recoil 
TA 


125 — 13 1.35 14 145. 


Mo, (GeV/c?) 


TOA 


FIG. 1: Scatter plots of Мол versus Мт! for (a) Ј/% and (b) 


(3686) data. The dashed lines denote the (1385)? signal region, 
and the solid lines denote the =° signal region. 
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FIG. 2: Distribution of Мол for (a) J/ and (b) (3686) data. 
The arrows denote the applied requirements, where the dashed arrows 
the ©(1385)° signal region and the solid arrows show the ЖО signal 
region. 


To determine signal yields, the mass of л°А is 
required to be within +34 MeV/c? for J/$ —> 


3:(1385)99:(1385)9, +10 MeV/c? for J/y — 2020, 
+35 MeV/c? for (3686) — X:(1385)95:(1385)9, and 
+11 MeV/c? for (3686) — =°=°, around the nom- 
inal mass of ©(1385)°/=°; the requirements are opti- 
mized by the FOM. For the (3686) decays, the re- 
quirements of |M?**?! — M j| > 0.005 GeV/c? and 
2991 — Mj j| > 0.015 GeV/c? are used to suppress 
the backgrounds (3686) — тт.//ф, where Mrecoil 
and Med are the recoil masses of any тїт” and 
толо combination if found, and M уу, is ће J/«» nom- 
inal mass according to the PDG [4]. 


4. Background study 


The data collected at CM energies of 3.08 GeV (30 
рь!) [14] and 3.65 GeV (44 рь!) [16] are used to es- 
timate the contributions from the continuum processes 
етет — X(1385)93(1385)" and =°=°. By applying 
the same event selection criteria, only a few events sur- 
vive and do not form any obvious peaking structures in 
the 4(1385)°/=° signal regions in the corresponding 
MD distributions. Taking into account the normal- 
ization of the luminosity and CM energy dependence of 
the cross section, the QED backgrounds are found to be 
negligible. 

The contamination from other background sources 
is analyzed using samples of MC simulated events of 
generic Y% decays that contain the same number of 
events as the data. After applying the same event se- 
lection, it is found that the peaking backgrounds for 
the  — X(1385)99:(1385)? mode mainly come from 
Vy — 2020, E-E*, (1385) £(1385)*, E(1530)E + 
c.c., and 70 12(1385)0, where the branching fractions 
for  — Е(1530)& + c.c. апа 19A3:(1385)? are taken 
from the isospin partner modes 7/0 — =(1530)= + 
c.c.[4] and т A3(1385) [13] based on the assumption 
of 12% rule. For ће 7/0 — =°=° mode, the peak- 
ing backgrounds are found to be from 7/0 — E-&*, 
уп. (39350, y£0 0), 3:051(1385)9, and 00, For the 
(3686) — =0=° mode, the peaking background is 
from (3686) — X939, and other backgrounds are 
found to be distributed smoothly in Mee mass spec- 
trum. 

The final states of baryon and anti-baryon decays 
both include a neutral pion with almost the same mo- 
menta. The 7° from the anti-baryon can be wrongly 
combined with the A in the ©(1385)°/=° reconstruc- 
tion. As a result, the wrong combination background 
(WCB) in the л°А mass spectrum is inevitable. This 
background is studied by the MC simulation. 


5. Results 


5.1. Branching fraction 


The signal yields for the decays v  — 
3:(1385)03:(1385)? and =°=° are extracted by per- 
forming an extended maximum likelihood fit to the 
MM spectrum. In the fit, the signal shape is repre- 
sented by the simulated MC shape convolved with a 
Gaussian function to take into account the mass res- 
olution difference between data and MC simulation. 
The peaking backgrounds and the wrong combination 
background are described by the individual shape taken 
from MC simulation, and the corresponding numbers 
of background events are fixed according to the individ- 
ual detection efficiencies and branching fractions [4]. 
The remaining backgrounds are found to be distributed 
smoothly in the AP spectrum and are therefore 
described by a second-order polynomial function. Fig- 
ure 3 shows the projection plots of Mr for the de- 
cays № — X(1385)93:(1385)? and =°=°, respectively. 

The branching fraction can be calculated by 


Bly > XX) = yeaa Bü EG 
where X stands for the X:(1385)? or =° baryon, є de- 
notes the detection efficiency obtained with the mea- 
sured o value, Nops is the number of observed signal 
events, B(X — «? A), B(A > pr) and B(T? — уту) 
are the branching fractions of X — Ал°, A — pr and 
по — yy taken from PDG [4], Ny 15 the total number 
of J'/1b or (3686) events [14, 16]. Table 1 summarizes 
the numbers of observed signal events, the correspond- 
ing efficiencies, and branching fractions for the various 
decays in this measurement with the statistic uncertainty 
only. 


5.2. Angular distribution 


The values of a for the four decay processes are de- 
termined by performing a least squares fit to the cos 0 
distribution in the range from —0.8 to 0.8, divided 
into 8 equidistant intervals for the decays (3686) — 
3:(1385)93:(1385)? and into 16 intervals for the other 
three decay modes. 

The signal yield in each соз bin is obtained with 
the aforementioned fit method. The distributions of the 
efficiency-corrected signal yields together with the fit 
curves are shown in Fig. 4. The o values obtained from 
the fits based on Eq. (1) are summarized in Table 1. 
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FIG. 3: (Color online). Recoil mass spectra of m°A for (a) J/y — %(1385)°4(1385)9, (b) J/b — EZ, (c) (3686) — 
©(1385)95(1385)°, and (d) 4/(3686) — =0=0. Dots with error bars indicate the data, the blue solid lines show the fit result, the red short-dashed 
lines are for signal, the red long-dashed ones are for the remaining background (Other-Bkg), and the green hatched ones are for wrong combination 
background (WCB) , the black hatched ones are for the peaking backgrounds. 


TABLE 1: The numbers of the observed events Nops, efficiencies є, о values, and branching fractions B for y — X:(1385)9 (1385)? and 2020, 
Only the statistical uncertainties are indicated. 


Channel Nobs e(96) a В(х10-%) 
Jb > У:(1385)05(1385)0 102762 + 852 13.32 + 0.04 —0.64+0.03 10.71 + 0.09 
Ј/р — B=? _ 134846 + 437 14.05 + 0.04 0.66 + 0.03 11.65 + 0.04 
(3686) > x:(1385)92:(1385)? 2214 + 149 13.13 + 0.03 0.59 + 0.25 0.69 + 0.05 
(3686) — =0=0 10839 + 123 14.10 + 0.04 0.65 + 0.09 2.73 + 0.03 

6. Systematic uncertainty range, wrong combination, the background shape, and 


the angular distributions. In addition, the uncertainties 


6.1. Branching fraction of the decay branching fractions of intermediate states 


Systematic uncertainties on the branching fractions and uncertainties of the total number of v events ar 
are mainly due to efficiency differences between data also accounted for in the systemate uncertainty, АП of 
and MC simulation. They are estimated by c ompar- the systematic uncertainties are discussed in detail be- 
ing the efficiencies of photon, 7°, A and E? recon- low. 
struction between the data and the MC simulation. Ad- 1. The uncertainty associated with photon detec- 
ditional sources of systematic uncertainties are the fit tion efficiency is 1.0% per photon, which is de- 


termined using the control sample 7/0 — pm. 
Hence, for  — 4(1385)°X(1385)°, the value 
2.096 is taken as the systematic uncertainty. 

. The systematic uncertainty due to the 1C kine- 
matic fit for the 7° reconstruction is estimated to 
be 1.0% with the control sample J/4» — pr. 


. The uncertainty related to the A reconstruction 
efficiency in (1385) decays is estimated using 
the control sample 7) — =~=+. Here, the A re- 
construction efficiency includes systematic uncer- 
tainties due to tracking, PID, and the vertex fit. A 
detailed description of this method can be found 
in Ref. [22]. 
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FIG. 4: Distributions of cos@ for (a) J/v  — 
3:(1385)0$(1385)9,(b) J/4 — &989, (c) (3686) > 
3:(1385)9 $:(1385)9, and (d) (3686) — =0=0. The dots 
with error bars indicate the efficiency corrected data, and the 
curves show the fit results. 


. The E? reconstruction efficiency, which includes 
the two photon efficiencies, 1? reconstruction ef- 
ficiency and the A reconstruction efficiency, is 
studied with the control sample //ф — =°=° 
via single and double tag methods. The selection 


criteria of the charged tracks, and the reconstruc- 
tion of А and E? candidates are exactly same as 
those described in Sec. 3. The =° reconstruction 
efficiency is defined as the ratio of the number of 
events from the double tag =°=° to that from the 
single tag. The difference in the E? reconstruc- 
tion efficiency between data and MC samples is 
taken as the systematic uncertainty. 


In the fits of the Мт! signal, the uncertainty 
due to the fitting range is estimated by varying the 
mass range by + 10 MeV/c? for two sides. The 
resulting differences of signal yields are taken as 
the systematic uncertainty. 


. The uncertainties due to the background shape 


arise from the polynomial function and the peak- 
ing shape. The former is estimated by the alter- 
native fits with a first or a third-order polynomial 
function. The latter is estimated by varying the 
number of normalized events by 1c. The larger 
difference is taken as the systematic uncertainty. 
The total uncertainty related to the background 
shape is obtained by adding all contributions in 
quadrature. 


. The systematic uncertainty due to the wrong com- 


bination background is estimated by comparing 
the signal yields between the fits with and with- 
out the corresponding component included in the 
fit. The differences of signal yields are taken as 
systematic uncertainties. 


. The uncertainty related with the detection effi- 


ciency due to the modeling of the angular dis- 
tribution of the baryon pairs, represented by the 
parameter a, is estimated by varying the mea- 
sured а values by Io in ће MC simulation. The 
changes in the detection efficiency are taken as a 
systematic uncertainty. 


The systematic uncertainties due to the branching 
fractions of the intermediate states, =°, X:(1385)? 
and A, are taken from the PDG [4]. They are 
1.9% for у — X:(1385)93:(1385)? and 0.8% for 
y — 5050. 


. The systematic uncertainties due to the total num- 


ber of J/w or (3686) events are determined 
with the inclusive hadronic v» decays. They are 
0.5% and 0.6% in [14, 16], respectively. 


The various systematic uncertainties on the branch- 


ing fraction measurements are summarized in Table 2. 
The total systematic uncertainty is obtained by sum- 
ming the individual contributions in quadrature. 


TABLE 2: Relative systematic uncertainties on the branching fraction measurements (in 96). 


Jh > (3686) — 
Source $(1385)0x(1385) =020 x(1385)X(1385)0 2020 
Photon efficiency 2.0 2.0 
по reconstruction 1.0 1.0 
А reconstruction 3.0 E 1.0 E 
=° reconstruction 2. 2.6 " 2.6 
Fit range 2.1 1.6 2.8 1.8 
Background shape 3.9 1.5 4.0 2:3 
Wrong combination 4.2 0.8 4.5 0.3 
Angular distribution 2.0 0.5 1.2 2.8 
Intermediate decay 1.9 0.8 1.9 0.8 
Total number of 4 0.5 0.5 0.6 0.6 
Total 73 3.7 7.4 4.9 


6.2. Angular distribution 


Various systematic uncertainties are considered in 
the measurement of the values of œ. These include the 
uncertainty of the signal yield in the different cos 0 in- 
tervals, the uncertainty of the cos 0 fit procedure, and the 
uncertainty related to the detection efficiency correction 
curve as function of the cos 0 bin. They are discussed in 
detail below. 


1. The signal yields in each cos @ interval are de- 
termined by the fit to the corresponding Ap 
distribution. The sources of the systematic un- 
certainty of the signal yield include the fit range, 
the background shape, MC resolution and wrong 
combination, where the MC resolution is fixed 
for the decay (3686) — X(1385)93:(1385)9 
only. To estimate the systematic uncertainty re- 
lated with fit range on Мт, we repeat the fit to 
the M. iy dm distribution by changing the fit range 
by +10 MeV/c?. Then, the o values are extracted 
by the fit with the changed signal yield, and the 
resulting differences to the nominal a values are 
taken as the systematic uncertainties. The uncer- 
tainties related to the background shape, MC res- 
olution and wrong combination backgrounds in 
the fit are evaluated with a method similar to the 
one described above. 


2. The systematic uncertainties related to the proce- 
dure of the fit on the cos distribution are esti- 
mated by re-fitting the cos distribution with a 
different binning and fit range. We divide cos 0 
into 8 intervals for у — 2020 and 16 intervals for 
V — %(1385)°S(1385)°. The changes of the a 


values are taken as systematic uncertainties. We 


also repeat the fit by changing the range to [—0.9, 
0.9] or [—0.7, 0.7] in cos 6, with the same bin size 
of the nominal fit. The largest differences of a 
value with respect to the nominal value are taken 
as a systematic uncertainties. 

3. Inthe analysis, the o values are obtained by fitting 
the cos distribution corrected by the detection 
efficiency. To estimate the systematic uncertainty 
related to the imperfect simulation of the detec- 
tion efficiency, the ratio of detection efficiencies 
as function of cos@ between data and MC sim- 
ulation is obtained based on the control sample 
J/w = 2°=° with a full event reconstruction. 
Then, the efficiency corrected соз distribution 
scaled by the ratios of detection efficiencies is re- 
fitted. The resulting differences in o are taken as 
the systematic uncertainty. 


АП the systematic uncertainties for the œa measure- 
ment are summarized in Table 3. The total systematic 
uncertainty is the quadratic sum of the individual values. 


7. Conclusion and discussion 


Using (1310.6 7.0) x 10° 7/ and (447.9+2.9) x 
10° (3686) events collected with the BESIII detec- 
tor at BEPCII, the branching fractions and the angular 
distributions for y — %(1385)°S(1385)° and =0=0 
are measured. A comparison of the branching frac- 
tions between our measurement and previous experi- 
ments (PDG average) is summarized in Table 4. The 
branching fractions for Y — Х(1385)0(1385)° are 
measured for the first time, and the branching fractions 
for ф — =°=° are measured with a good agreement 


TABLE 3: Relative systematic uncertainties on the a value measurements (in 96). 


Jf > (3686) > 
Source $(1385)0x(1385) =0=7 (1385) X(1385)0 2020 
aen fitting range 7.8 3.0 15.3 77 
Background shape 3.2 3.0 20.0 4.6 
MC resolution Т se 16.9 as 
Wrong combination 4.7 1.5 5.1 15.0 
cos 0 interval 7.8 3.5 22.0 10.4 
cos @ fitting range 7.8 3.0 15.6 3.5 
Efficiency correction 4.7 3.0 9.0 3.0 
Total 15.4 7.1 41.8 20.8 


and a much higher precision than the previous results. 
The measured a values are also compared with the pre- 
dictions of the theoretical models from Refs. [6, 7]. As 
indicated in Table 5, some of our results disagree sig- 
nificantly with the theoretical predictions, which may 
imply that the naive prediction of QCD suffers from 
the approximation that higher-order corrections are not 
taken into account. As calculated in Ref. [9], the sign 
for parameter a in i — X93? mode could be negative 
if re-scattering effects in the final states are taken into 
account. However, our results show that o for J/4b is 
negative, and is different to the other decay processes in 
this measurement, which is hard to explain within the 
existing models. We, therefore, believe that it is of ut- 
most importance to improve the theoretical models to 
shed further light on the origin of these discrepancies. 


To test the “12% rule", the ratios of the 


: : B(1 (3686) —5X:(1385)9 3:(1385)?) 
branching fractions C BJS (1385)95(1385)0) - nd 
B(v(3686) 5 £92? 


—BUJG EU. are calculated to be (6.44 + 0.47 
0.64)% and (23.43 + 0.26 + 1.09)%, respectively, tak- 
ing into account the cancelation of the common system- 
atic uncertainties. The ratios are not in agreement with 
12%, especially for the =°=° final state. 

To test isospin symmetry, the ratios of the branch- 
ing fractions listed in Table 6 are also calculated based 
on the measurements between the neutral mode and the 
corresponding charged modes [13] taking into account 
the cancelation of the common systematic uncertainties. 
АП ratios are within 1c of the expectation of isospin 
symmetry. 
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